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Abstract 

Factorization in gauge theories holds at the amplitude or amplitude-squared level 
for states of given soft or collinear momenta. When performing phase-space integrals 
over such states, one would generally like to avoid putting in explicit cuts to separate 
soft from collinear momenta. Removing these cuts induces an overcounting of the soft- 
collinear region and adds new infrared-ultraviolet divergences in the collinear region. In 
this paper, we first present a regulator-independent subtraction algorithm for removing 
soft-collinear overlap at the amplitude level which may be useful in pertubative QCD. 
We then discuss how both the soft-collinear and infrared-ultraviolet overlap can be 
undone for certain observables in a way which respects factorization. Our discussion 
clarifies some of the subtleties in phase-space subtractions and includes a proof of the 
infrared finiteness of a suitably subtracted jet function. These results complete the 
connection between factorized QCD and Soft-Collinear Effective Theory. 
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1 Introduction 


Factorization is at the heart of our ability to use perturbative quantum chromodynamics 
(QCD) to make theoretical predictions for scattering processes at high-energy particle col¬ 
liders. It is extremely fortuitous that accurate particle distributions can be computed by 
convolving universal parton distribution and hadronization models with perturbative calcu¬ 
lations of jet formation. While factorization at the non-perturbative level is hard to establish, 
factorization relevant to the structure and substructure of jets can be understood within per¬ 
turbation theory. In particular, the radiation patterns in perturbative QCD factorize into 
hard, collinear and soft contributions. Moreover, subtleties in perturbative factorization (for 
example, related to non-global logarithms [1-7]) are a limiting factor in many ultra-precise 
jet-substructure calculations. Thus, there has recently been renewed interest in studying 
factorization, particularly in the context of Soft-Collinear Effective Theory (SCET). 

A concise formulation of factorization in QCD was proposed and proven in [8] and [9], 
hereafter referred to as [FS1] and [FS2] respectively. These papers build upon decades of 
insight [10-16]. Up to color factors, the formula from [FS2] reads: 


(x 1 ---x N] x a \o | 0 ) 


r-~j 


CoiSii) 


(W | ipW\ |0) 

<0|uM|0> 


(X N \W^\0) 
(0| w' n Yn |0> 


(Wlif-.-Uvlo) 


( 1 ) 


In this expression, the state {X\ ■ ■ ■ X N \ X s \ has soft particles, in (X s [, and particles collinear 
to various specified directions, in (X, \. The left-hand side is a matrix element in QCD of an 
operator like O — ip ■ ■ ■ ip in this state. The right hand side is a factorized product of matrix 
elements, each of which involves only one collinear sector or the soft sector. The symbol = 
indicates that the two sides are identical at leading power. More precisely, if one were to 
compute some infrared-safe observable dominated by soft or collinear radiation, such as the 
sum of the jet masses r = ^ °f the terms in 'jpp that are dominant as r —)■ 0 will 

be identical on both sides. More details can be found in Section 2 below and in [FS1] and 
[FS2], 

The formula in Eq. (1) presupposes that the external momenta are designated as soft or 
collinear. If a particular momentum can be classified as soft or collinear, then the factorized 
formula will hold whether it is put in (X s [ or in the appropriate (X t [. For example, we can 
place all the soft-collinear momenta in the soft sector by designating any particle with energy 
less than some A as soft, and then draw cones of size R around each of the hard directions 
to distribute particles in the collinear sectors. With such hard cutoffs, one can then square 
the matrix elements on the right-hand side of Eq. (1) and perform the phase-space integrals 
over the appropriate measurement function to get a differential distribution. The result will 
agree at leading power with the the distribution computed using the left-hand side of Eq. (1) 
in the limit R —>• 0 and A —>■ 0. 

There are two problems with the hard-cutoff prescription for resolving the soft-collinear 
ambiguity. The first is practical: introducing an extra scale makes the relevant calculations 
nearly impossible. Moreover, the cutoff dependence may not exactly cancel in the factorized 
expression and therefore one must either take R —> 0 and A —> 0 after the calculation or 
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live with power corrections in these cutoffs. The second is conceptual: the cutoffs violate 
factorization in the following sense. There will in general be leading-power dependence on the 
cutoff in the soft and collincar sectors separately (terms like f In R, for example) which only 
cancel when the sectors are combined. Thus the two sectors are not completely separated. 

It would be great if we could simply perform phase-space integrals over each sector 
separately including all momenta. This is not as crazy as it sounds. We know that including 
very energetic virtual momenta in the soft or collincar sectors causes no problem, since 
the modification can always be compensated for in the matching coefficient (Co(*S'p) in 
Eq. (1)). Indeed, effective theories always have different ultraviolet (UV) structure from the 
full theories to which they are matched. For example, in SCET, there are \ UV poles at 
1-loop in dimensional regularization, while in full QCD, one only ever has \ poles. In fact, 
these double poles allow for the resummation of Sudakov double logarithms in SCET using 
the renormalization group. We also know that one does not have to distinguish soft from 
collinear momenta in loops when using Eq. (1): the overcounting is compensated for by the 
vacuum matrix elements in the denominator of this equation. Thus, we have good reason to 
believe that subtractions similar to the denominator factors in Eq. (1) can be added to this 
formula to allow for unrestricted phase-space integrals. 

Removing the overcounting of soft and collinear momenta has been addressed in the 
traditional approach to factorization, for certain observables [17-19]. There, the soft limit 
of collinear momenta is compensated for with eikonal jet functions [20]. In SCET, the 
overcounting can be formally avoided by not including the zero-momentum bin in any of 
the collincar sectors [21], This exclusion translates into a subtraction diagram-by-diagram. 
This zero-bin subtraction is necessary in SCET because the same soft-collincar momentum 
region in QCD is represented by multiple fields in the effective theory (similar overcounting 
is present in other effective theories, such as NRQCD). In [22-24] the two prescriptions were 
shown to be equivalent. Alternatively, in the method-of-regions approach to SCET [25-28] 
the overcounting is sidestepped through careful consideration of the analytic properties of 
the contributions from different sectors. We briefly review these approaches and contrast 
them with our approach in Section 3.3. 

The formulation of factorization in [FS1] and [FS2] and Eq. (1) is intermediate between 
traditional QCD and SCET. It provides a precise formulation of factorization purely in 
terms the holds in full QCD, but has a factorized form with a natural effective held theory 
interpretation. It is based on the observation of Freedman and Luke [16] that the unwieldy 
Feynman rules of SCET can be avoided and the effective Lagrangian taken simply as the 
direct sum of N + 1 copies of the QCD Lagrangian, corresponding to N collincar sectors and 
a soft sector. The formulation in [FS1] and [FS2] can be thought of as a generalization of 
the Freedman-Luke proposal, equivalent but simpler at leading power, and that addresses 
the soft-collincar overlap of virtual momenta. In this paper, we extend the formulation so 
that phase-space integrations can be done without explicit cutoffs on the momenta of various 
sectors. 

There are two main results in this paper. First, in Section 2, we show how the specification 
of which sector a gluon belongs to can be removed at the amplitude level. More precisely, 
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suppose we have an amplitude Ai(pi, ■ ■ ■ ,p n , Qi • • • q m ) with n hard momenta and m other 
momenta in QCD. We show how an approximation to JA which we call A4 su b can be derived 
with the property that when any of the g* become soft or collincar to any of the pi, A4 su b 
agrees with JA at leading power. That is, one does not have to specify which sector the qt 
belong to - the matrix element is correct no matter what. While A4 su b is not a factorized 
product of matrix elements, it is the sum of factorized products of matrix elements of fields 
and Wilson lines. Each term in this product is simpler than full QCD. Thus such a subtracted 
matrix clement may be integrable analytically and therefore provide a useful basis for a 
subtraction scheme in fixed-order QCD. 

The second result, in Section 3, is a derivation of how at the amplitude-squared level 
factorization can be preserved and phase-space cutoffs removed for certain inclusive event 
shapes. Although the result of this section agrees with the eikonal-jet function subtraction 
method of traditional QCD (which is itself equivalent to SCET), we believe our derivation 
elucidates some subtleties and makes the procedure more systematic. In addition, we present 
explicit 1-loop formulas for various relevant soft and jet functions, with and without cutoffs 
and with different regulators. These formulas demonstrate which objects are infrared safe, 
cutoff-dependent, and well-defined. Some calculational details are relegated to the appendix. 
Section 3.3 contrasts our approach with previous approaches. We conclude in Section 4. 


2 Factorization at the amplitude level 


We begin by quickly reviewing the notation and main results of [FS1] and [FS2], These 
papers showed that factorization holds for massless particles whose momenta are either soft 
or collinear to one of N directions rJ-. States with particles of these momenta are written 
as (Ad ■ ■ ■ Xn] X s \. The hard scale (such as the center-of-mass energy) is denoted as Q and 
scaling parameters A j are defined each collinear sector and A s for the soft sector. Momenta 
in each collinear sector scale as 


(Xj | (■ ■ ■, qj , ■ ■ ■ I ^ 

and momenta in the soft sector scale like 


Q K • ( lr "J ■ 'lr % ) ~ (A). 1, Xj) 


(As | (• ■ ■, k s ,... | 


q j ’ fcsi Tij ' ksj ~ (A s , A s , As), Vj 


( 2 ) 

(3) 


For simplicity, assume the scattering process under consideration is the decay of a heavy 
particle mediated by an operator O in QED (to avoid cumbersome color indices of QCD). 
Then, the factorization formula takes the form of Eq. (1): 


(X 1 ---X N] X s \O\0) 




CoiSij) 


(XiIVTFQO) 

(o|rM|o) 


(X N \W^\0) 
(0| W ] n Y n |0) 


(As|y 1 t ---y iV |o) 


Here Co(Sij) is a finite function of the large products of the net momentum in each jet, 
Sij = P{ ■ Pj ; it does not depend on the small power-counting parameters, A j of A s . The 
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( 4 ) 


Wilson lines, Wj and Yj. are defined in QCD as follows: 



exp 


W 


dsrtj ■ A(x + srij ) e" 


and 



exp 


i9 



dstj ■ A(x + stj ) e 



(5) 


where t'- are some lightlike directions assumed not collincar to their associated n 1 - . The P{} 
denotes path ordering; in QED the path ordering is trivial and the electromagnetic charge 
is e = —g. Eq. (1) is an equality at leading power in all of Xj and X s separately. For many 
applications, such as for thrust, one takes A^ = A s for all j; in the SCET literature, this power 
counting is referred to as SCETi [29]. For recoil sensitive observables like jet broadening, one 
takes Xj = X s as in SCETn [30]. The factorization in Eq. (1) holds for any relative scaling. 

The important physics contained in Eq. (1) is that each factor on the right-hand side 
represents a different factorized sector: the Wilson coefficient, Co(5V,), represents all of the 
hard physics and must be IR-insensitive. Each collinear sector is represented by the ratio 
(Xj Wpb |0) / (0| WjYj |0) and contains only rq-collinear IR divergences. Finally, the soft 
sector is fully described by the matrix element, (AA.| Y-± ■ ■ ■ Y N |0), which contains all of the 
soft divergences of the full amplitude on the left-hand side of Eq. (1). 

One attractive feature of Eq. (1) is that each matrix element is constructed out of full- 
theory operators and evaluated using the full-theory Lagrangian; there are no additional 
subtractions/prescriptions needed, just simple QCD/QED Feynman rules. Moreover, the 
power counting is a consequence only of the scaling of the external momenta in the states 
(.Xi ■ ■ ■ X N : X s \ O |0). An obvious fact with important repercussions is that Eq. (1) is not 
valid when any of the momenta in a given sector does not obey the scaling that is associated 
with that sector. Consequently, one cannot, for example, integrate over the entire phase 
space of one of the external momenta in Eq. (1) because it would enter the scaling regime 
of other sectors. 

Therefore, when calculating cross sections by squaring Eq. (1) one can either integrate 
over the phase space dlixj with cutoffs in the integrals restricting each integral to be within 
the collincar region, or one can try to extend the integrations to the entire phase space and 
perform a subtraction that gets rid of the errors that we introduced by extending dlix to the 
entire phase space. Introducing cutoffs to integrals is incredibly tedious and produces new 
scales in the effective theory that obscure factorization (as shown explicitly in Section 3). 
The subtraction procedure is the only reasonable way forward. We next discuss subtractions 
at the amplitude level, and discuss subtractions at the cross section level in Section 3. 


2.1 Example subtractions 

Consider the case of a qqg final state, with quark momenta p V and p% in different direc¬ 
tions and the gluon momentum q A Suppose we want to integrate over the gluon momenta 
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( 6 ) 


inclusively. We can do so using Eq. (1) if when q || p 1 we use 


M l {y 1 ,p 2 ,q) 


wbm vW lm 

(0| y!w, |o> (o| w\y 2 |o> 1 2 


if q || p 2 , we use 


M 2 (pi,p2,q) 


(pil^W iJO) (p^q\Whm 

(o|rM|o) (o| w\y 2 |0) 


(o|nV 2 |o) , 


and if q is soft, we use 


M s (pi,p 2 ,q ) 


(o|yM|o) (o| w}y 2 \o) 


(7) 


( 8 ) 


However we split up the integration regions (say with a soft energy cutoff A and cone radius 
R) the dependence on the split (on A and R ) will drop out at leading power when all three 
contributions are added. Nevertheless, it would be nice to have an expression that we could 
simply integrate over q without ever introducing A and R in the first place. 

To proceed, we first examine the consequences of soft-collinear factorization for the op¬ 
erator (D^ w = ijjWi (rather than a local QCD operator like -0-0). The all-orders proof of 
factorization in [FS2] applies to 0^ w . In particular, if we have a state with momenta p\ ■ ■ ■ p n 
all of which are collinear to each other as well as momenta q\ ■ ■ ■ q rn all of which are soft, 
then 

(Pi---Pn,qi---qm\^W 1 \0) = Cj W (Pi---Pn\ W 1 10) ^ ^ (9) 

for some C^ w . To determine C$ w , we note that C^ w does not depend on how the momentum 
in the collinear and soft sectors are distributed; this equation holds for any n > 0 and any 
m > 0. In particular, if we take m — 0 then the two sides are identical (and agree at leading 
power) if and only if C^ w = 1. Thus we must have Cj, w = 1 for any states. 

As a special case, Eq. (9) implies that for one collinear and one soft momentum 


(PHgl^WilO) (PilVAWlO) (q\ Y^W\ |0) 

(o| fM |0) (0| Yjw 1 10) (0| fM |0) 

Similarly, applying the general factorization formula to O = F/F 2 , we get 


(9|^iV 2 |0) 


<? II pi 


(g|FM|0) 

(oifMIo) 


x (oif/f^o) 


( 10 ) 


( 11 ) 


In this case one can see that the Wilson coefficient is 1 to all orders by using the proof in 
[FS2] that the factorization theorem is independent of the collinear Wilson-line direction, f 1; 
and then choosing = n 2 , so that W\ = F 2 . 

With these results, we can now analyze the following all-loop-order subtracted matrix 
element: 
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M suh (pi,p 2 ,q) 


(pr q\ jiWi |o) _ (pi \ iJjWi |o) (g|yM| 0 ) l (palwftHo) yty 
(o|yM|o) (o|yM|o)<o|yM|o) J ( o \ w } y 2 \ o ) u 1 21 


(pi| TpWi |0) / ( P 2-,q\^W 2 \0) 

(0| YfWi |o) \ (o\y]w 2 |0) 


(p2\^w 2 | o ) (g\Y^w 2 \o) 

(o|y 2 V 2 |o) (o|y 2 V 2 |o) 


(o| yJy 2 |0) 


+ 


W«Wf8 w v,| 0) 

<o| y/w, |o> (o| w 2 f y 2 |o> 


(12) 


If we take q soft, then neither of the first two lines contribute by Eq. (10), and the result 
is given by the third line which is the correct leading power matrix element Ai s . When 
q || pi, then neither term in the second line is IR sensitive and the subtraction term in 
the first line (which is IR-sensitive) is canceled by collinear limit of the third line, using 
Eq. (11). Thus, only the Erst term on the first line contributes at leading power in this limit, 
in agreement with Ai\. The analogous argument works for the q || p 2 limit. We conclude 
that M suh (pi,p 2 , q) agrees with full QCD at leading power for any q. Thus, we can integrate 
-M su b over phase space without splitting the soft and collinear sectors. 

To be explicit, we can evaluate Eq. (12) in perturbation theory. At tree-level, 


M(pi,p 2 ,q ) = u(pi) 



Ml + 0 + gtl ■ RO 
2pi- q ti ■ q J 

f-2 + + ~Qt 2 ' £g \ _ 

2 p 2 - q t 2 - q ) 


f -gni • e q gti • e q 
\ 7ii- q ti ■ q 


v{P2 ) 


f ~gt 2 • e q gn 2 • e q 
\ t 2 - q n 2 -q 



+ u(pi) 


f -gni • e q 
V rii ■ q 


+ gn 2 • 
n 2 ■ q J 


v(p 2 ) (13) 


where each term in round brackets corresponds to one of the matrix elements containing the 
gluon, and thereby each satisfies the Ward identity separately. From the explicit expression 
in Eq. (13) it is easy to check that each soft and collinear limit works out exactly as stated 
in the paragraph after Eq. (12). It can also be seen that the tj dependent terms cancel out 
completely as do the soft terms containing rij at this order, leaving: 


M SVLh (pi,p 2l q) *= u(pi) 


-9t q ty i + j) 

2pi ■ q 


+ 


+ 


2 p 2 ■ q 


v{p 2 ) =° M(pi,p 2 ,q) 


(14) 


So the the full matrix element of QED is reproduced exactly in this case. Of course, for 
more complex calculations we expect A4 to only reproduce the full-theory matrix element 
at leading power, rather than be exactly equal to it. 
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2.2 General amplitude-level subtraction 

The generalization of Eq. (12) for arbitrary collincar and soft sectors is 


(X 1 ---X N -X s -q\O\0) 


(giljbdjO) 

(0|rM|0) 


(X N \ wy |0) 

(o| W' N Y N |0) 


{x s , q\ Yi ■ ■ ■ y n |o) 


f (X u q\W^\0) \ (X N \W^\0) , 

h <0|VM|0> l (0| W}Yi |0) /soft (0| WffYff |0) ( sl 1 


( 15 ) 


where the {jsoft notation means the operator matrix element corresponding to having sub- 

J sub 

tracted the q —>■ soft limit. To be explicit, we can use the notation S(q) as in [FS2] for the 
leading order contribution in the q —> soft limit. Then 


1 (Xi, y| Wjip 0) ' 

1 _(Xi,q\Wty\0) 


l (01 W}Y, |0) . 

[ A (0| W}Y, |0) 

V (o| w/y |o) ) 


_ (XiAwht iq) _ (vrgtyo) Mithjjg) (16) 

(o| w]y, |o) (o| w}y, |o) 2 

This subtracted quantity is exactly the same as what was used in Eq. (12) and vanishes at 
leading power in the q —>- soft limit by Eq. (10). Eq. (15) is a sum of factorized expressions 
which agrees at leading power with full QCD in any soft or collinear limit of q. 

To generalize to multiple gluons or quarks with momenta qi, the analogous formula is 
easiest to define recursively. For example, adding a second gluon to Eq. (15), we can either 
place it in the soft matrix element, or in a collinear matrix element. If it is in the collincar 
matrix element, we must subtract off the soft limit. Thus we get a sum of terms: 


(Xi • • • X N ] X s ; qi, q 2 \ O |0) = IR 


(XiIVTfdlO) 

(o|yM|o> 


(X N \ |0) 
(0| W ] n Y n |0) 


(X s , qi ,q 2 \Y?---Y N \0) 


tt l (0| W-Yi |0) J«ft 


(x s , q 2 \ Yj ■ ■ ■ yjv |o) 


N 



i,j= 1 


A f (X,, q 2 \ W'l'ijj |0) | 

tt l <o| w}y, |o) IS, 


<x„ 9 ,|y 1 , ...y K |o) 


l <0| W}Y, |0) 


<o wJy,- |0) 


I soft 
y sub 


(17) 


where the • • • represent the other collinear matrix elements which do not contain any g’s. 
In the last line, when i = j the soft subtraction must be done iteratively to ensure that the 
subtraction mitigates the soft enhancement in any order of limits of q± and q 2 going soft. 


















That is, 


{X,, qi ,qi\W]i,\0) 


<0| WjVj |0) 


91,92 = 
soft 
sub 


( 0 | W}Yj | 0 ) 


wim 

(0| w]y s |0) 


5 ( 91 , 92 ) 


(X J:gi ,q 2 \W^ |0) 

(0| WjYj |0) 


{X^quqtlWtylO) 


WjYj | 0 ) / 5(91, 92)/ 5 ( 9l ) 


'(X,, gi , g 2 | WU |0> / (X,, gi , g 2 | WU |0) 


(0|^|0) 


( 0 | | 0 ) / 5(91,92)/ 


( 18 ) 


where S(qi,q 2 ) means taking the leading-power expression in the qi,q 2 —> soft limit simul¬ 
taneously and, therefore, does not drop q\ with respect to q 2 or vice-versa. Note that, as 
always, we can write the soft limits in terms of amplitudes with Wilson lines using the 
factorization theorem of Eq. (1). For example, 


/ (y, g i,i 2 Hy>|o> \ = (V|ty>|o) fa.^ln/yio) 

V <0| WjYj |o) T„„, 2) (0| W}Yj |o> <o| wJy,- |0) 

where we know that the Wilson coefficient will always be 1 to all orders by the argument 
given after Eq. (9). 

That these subtractions will always work follows using the arguments of [FS2], In par¬ 
ticular, the “coloring algorithm” in Section 6 of that paper is exactly the recursive soft 
subtraction procedure indicated by Eqs. (12), (15)—(18) . As with the algorithm in [FS2], 
the soft limit of any subset of the g’s in the {} 5 soft '* matrix elements are power suppressed, 
and they should correspondingly be colored blue. With this knowledge, it is easy to check 
that Eq. (17) agrees in the IR: when q±, q 2 —> soft, all of the {/soft matrix elements are power 
suppressed and only the top line survives, which gives the correct answer. When q\ || pj and 
q 2 —» soft, say, the bottom two lines are power suppressed and the 5(gi)-subtracted term 
cancels with the top line, leaving only the one term that matches the full-factorized formula 
in this limit. Similarly, all other limits can be simply checked. The pattern of subtractions 
with more than two gluons follows exactly as with the coloring algorithm stated in generality 
in [FS2], 

The procedure outlined in this section produces amplitudes which can be computed as 
a sum of factorized terms. These amplitudes, which are a new result, reproduce all of the 
leading-power IR-sensitive limits of the full-QCD amplitudes, at all-loop order. Each fac¬ 
tor in each term in the sum involves matrix elements of fields and Wilson lines that are 
universal and simpler than the factors in the full QCD amplitude. Given these properties, 
an interesting application of the matrix elements derived in this section might be towards 
subtraction procedures for QCD calculations at NNLO or beyond. One application of sub¬ 
traction methods is to split an amplitude into a universal IR-sensitive piece that is simple 
enough to integrate analytically and a piece that is IR-hnite which could be integrated nu¬ 
merically [31-35]. The amplitudes presented in this section could be a candidate for such a 
procedure at any order in perturbation theory and for any number of external particles. 
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3 Factorization for distributions 


Despite having many strengths, amplitudes as in Eq. (15), are no longer factorized: they 
cannot be written as a single product of terms with the same external states (in this case 
the collincar sectors and the soft sector are tangled). When the amplitudes are squared, the 
interference effects between various terms in the sum contribute at leading power, so they 
must all be included. Thus, while one can integrate over the momenta qj without overcount¬ 
ing the infrared-sensitive region, the separation between soft and collinear contributions is 
no longer manifest. Moreover, it is not clear how the large logarithms associated with the 
leading-power 1R sensitivity can be resummed using such amplitudes. 

Fortunately, for certain observables, one can perform subtractions differently so that 
factorization is preserved at the cross-section level. In this section, we discuss a class of 
factorizing observables. Namely, we discuss observables whose measurement function, that is, 
the mapping from the final-state momenta to the observable, is linear in the soft and collincar 
momenta. These observables include many e + e _ event shapes, such as thrust [?, 36-39], 
angularities [20,40,41] heavy jet mass [42], the C parameter [43-45] and jet broadening 
[46-50]. Many hadron collider observables are also in this class [51], such Drcll-Yan near 
threshold [52], deep inelastic scattering as x —> 1 [53], direct photon production [54,55], W/Z 
+ jet [56-58], jet mass [59,60], or tt production near the hadronic threshold [61,62] as well 
as lV-(sub)jettiness [63-65]. 

Factorization at the cross-section level for observables in this class has been understood 
already by traditional QCD and by Soft-Collinear Effective Theory (see above references). 
The overcounting of soft and collinear integration regions is also well-understood in both 
approaches, and the two approaches have already been shown to be equivalent [22-24], 
Unfortunately, it is challenging to extract from the literature which aspects of the removal of 
overcounting have been rigorously proven (in either approach) and which aspects are simply 
assumed. Moreover, the overcounting in phase-space integrals has not been addressed at all 
in the effective field theory formulation with full-theory fields [16], [FS1] , [FS2] . The goal 
of this section is to give a self-contained proof that the overcounting induced by removing 
phase-space cutoffs can be completely compensated for. We thereby demonstrate a form of 
factorization that holds exactly at leading power at the cross-section level with no phase 
space cutoffs. 


3.1 Factorization for thrust 


For concreteness and simplicity, we begin our discussion with thrust, the paradigmatic ob¬ 
servable whose distribution factorizes. Thrust, T, is defined as [?] 


_ Ej I Pj ' n 

Sj \Pj\ 


( 20 ) 
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Figure 1: Mutually exclusive sectioning of phase space into j -collinear, soft and hard mo¬ 
mentum labelled by Rj, A R and //. respectively. In explicit formulas in this paper, R is 
treated as a rapidity variable: R = tan 2 |, with 6 the opening angle of the cone. 


where n is the thrust axis, defined to maximize T. The region where factorization holds is 
where r = 1 — T -Cl. Then 


T - ( 21 ) 

j * * j 

where Q is the center of mass energy and fl T (p) is the measurement function for thrust: 

Q T (p) = p~9(p + — p~) + p + 6(p~ — p + ) (22) 

where p + = n ■ p and p~ = n ■ p. 

Note that r has the property that it is linear in the momenta: each particle momentum 
contributes additively to thrust, independent of the other momenta in the final state (X\. In 
particular, if we decompose (A"| into soft, collinear and hard momenta, then we can compute 
the contribution to thrust from each sector separately and just add the results. In other 
words, linearity implies 



where p s is the sum of Vt T {k) over the soft momenta, p\ and P 2 the sum over collinear momenta 
in each direction and ph the sum over the remaining momenta. Writing the argument of the 
5-function as a sum lets us turn products of matrix elements into convolutions. 

To be concrete let us place the momenta into sectors using hard cuts: we draw cones 
of angular size R around the n and — n collinear directions and a ball of size A around the 
origin; anything in the cones but not the ball is collinear, (X,j, anything in the ball but not 
the cones is soft, (X s \. For later convenience, we include the soft-collinear radiation, which 
is in both the ball and the a cone, in the collinear sector (we could equally well have put it 
in the soft sector). Anything not in the cone or ball is called hard, (Xh |. This breakdown 
of phase space is shown in Fig 1. 
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Consider the thrust distribution in full QCD mediated by the operator O = That 

is, 

% = Y.J ^ x \{X\O\0)\H(t - ^n T (p x )) (24) 

where the sum is over all possible final states (X\ and the normalization and momentum- 
conserving 5-function are left implicit. 

When r is small, only states of the form (X = (Xj (X 2 | (X s | contribute at leading power 
in r. With the hard phase-space cuts in place, the factorization formula at the amplitude 
level, Eq. (1) along with Eq. (23), immediately generates a factorization formula for the 
thrust distribution: 

^ = H x S AR ® J Rl ® J R2 (25) 

dr 

Here H = \C\ 2 refers to the hard function (the square of the Wilson coefficient in the 
factorization formula). J Rl and J R 2 are jet functions with restricted phase-space integrals: 


" ,Rl M = E f 

A', J 


(WlVdWlO) 

(0|rM|0> 


S ( T ~Px 1 ) 


(26) 


with Xi the set of states all of whose momenta are within an angular distance R of the n 11 
direction. Note that we have modified the measurement function from f2 r (p) in Eq. (22) to 
simply p + . This is allowed since all the momenta in the cone necessarily have p~ > p + so 
the step functions in Eq. (22) can be evaluated explicitly. Analogously, the jet function, J R '\ 
will have the measurement function replaced by p~. Lastly, S AR is the phase-space restricted 
soft function 


S AR (r) = £ 

X s 


dRx a 



(27) 


Here the states have momenta which are not collinear, that is, they are an angular distance 
greater than R from all jets, and they have energy less than A. 

The equivalence in Eq. (25) holds at leading power in r only if R and A are small enough 
so that the collinear radiation is collinear and the soft radiation is soft. More precisely, it 
holds at leading power in R and A, meaning that the two sides my differ by terms of order R 
or order A which vanish as R —> 0 and A —y 0. Since the operators entering the soft and jet 
function are different, we do not expect the R dependence to cancel exactly between them; 
the factorization theorem only guarantees that it vanish at leading power. 

We suspect it may be pedagogically useful to examine explicit expressions for S AR and 
J R . To distinguish UV divergences from IR divergences we include an off-shellness regulator 
oj for the IR (see Eq. (50) in the Appendix) and analytically continue to d — 4—2e dimensions 
for the UV. Some intermediate steps are relegated to Section 6.2 in the Appendix. For the 
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unrenormalized jet function at finite R, we find 


J Rj (t) 2* S(t) + 


Oi s C F 


2tt \Q 


/j‘\ I , . / 2 3 7 7T 2 




6 


+ S(t) ( — - lnu; — 2 In a; In R + 2 In 2 a; + O(R) 


2 Ini? 


lnr 

r 


(28) 


This expression includes both the real and virtual contributions to {X\ \ i[)W\ |0) and the 
purely virtual contributions to (0| Y{W\ |0) in Eq. (26). Note that it has \ UV poles, which 
come from the virtual graphs. It also has an overlapping UV-IR singularity (the - In a; term 
on the second line). This singularity, which cannot be removed through local counterterms, 
comes from loops involving the Wilson lines which go to infinite energy collincar to one of 
the Wilson line directions. The log R dependence in Eq. (28) comes from the soft-collinear 
region of the restricted phase-space integral. Indeed, it cannot come from the collinear-but- 
not soft region, since at arbitrarily small r, the radiation is forced arbitrarily close to the jet 
axis and must be a finite distance from the cone boundary. In the soft-collinear region, the 
radiation can be soft but an angular distance R from the axis, so there can be R dependence 
at leading power in r. That the In R dependence comes from only the soft-collinear region 
is to be expected if it is to be completely canceled by the soft function. 

The soft function outside the cones with finite A is 


S ar (t ) = 5(t)+C f — ( ) { 5(t) ( X— + -\nu + 2hiu\nR-2ln z 0J + O(R) 


& S ( 
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£ Q 2 


2 tQ' 
- In — 
r A 




(29) 


This function also has an incomplete cancellation between the real and virtual contributions. 
In particular, the virtual includes the soft-collinear region which is excluded from the real 
emission. Note that the A dependence is entirely subleading power in r: for small r, the 9 
function in the third line in Eq. (29) vanishes and the other 9 function evaluates to unity. 
The i? dependence is not subleading power as r —> 0. There are also O(R) terms not shown 
here but written out in Eq. (82). We will come back to the cancellation of the i? dependence 
among the two jet functions and the soft function shortly. 
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Convolving Eq. (28) for each jet with Eq. (29) we get 


S AR ® J Rl ® J R -‘ 


+ (30) 

Note that the A dependence has dropped out completely, and the R dependence which is 
singular as R —* 0 has also dropped out. It is not hard to verify that this result agrees 
with the full-theory result for thrust at leading power, up to the coefficient of 5{r) which is 
corrected by the hard function. 

While the factorization formula for thrust in Eq. (25) works, it has numerous flaws. On 
the practical side, it is difficult to use because of the phase-space cuts. On the conceptual 
side, the cuts introduce additional scales into the soft and jet functions which frustrate 
factorization and resummation. The most serious flaw, however, is that the jet and soft 
functions are not individually infrared safe: they each have infrared divergences which cancel 
only when combined, as we saw with the explicit example above. These divergences come 
from an incomplete cancellation between the real-emission graphs, which have phase-space 
restrictions, and the virtual graphs, which do not. We could attempt to put phase-space 
cuts on the virtual graphs as well. However, it is more logical to try to remove the phase- 
space cuts from the real-emission contributions to the jet and soft functions, since this would 
simplify their calculation and removes the spurious scales. 

First we remove A. This is quite simple. Only the soft function depends on A. By our 
definition, S AR in Eq. (25) only integrates over the soft-but-not-collinear region of phase 
space. The phase-space region outside of the cones but with energy A < E < oo does not 
contribute at all at leading power in r. So we can simply define a new soft function by 
including also this region: 

S Ti 9* S AR (31) 

where S R = 5'0= 00 ) /? has no cutoff on energy in the soft function. This equivalence can be 
verified at order a s in Eq. (85), where the entire A dependence is subleading power in r, as 
observed above. Explicitly, 


S « T > 


6(t) [1 + 1 + 1-^ + 0(R)) - 5 


2e 


12 
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S R (r) =5{t) + C f ^ 
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T 

+ J 


(32) 


Note that S R is identical to the coefficient of 9 (A — ^) in Eq. (85). One might have imagined 
that taking A —>■ oo would introduce new UV poles. However, radiation in R, say in the 
right hemisphere, at a given r must have k + = Qt and k~ < j^k + = so the energy 
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E = \{k + + k ) < ^-(1 + jj) of all radiation contributing to S R at fixed r is in fact bounded 
from above so there are no new UV divergences. Thus, we now have 

— = H x S R ® J Rl <g) J R 2 (33) 

ar 

with no A dependence on either side. Keep in mind that this equivalence is still valid only 
as R —> 0: there are power corrections in R on the right-hand side. 

Removing the R dependence is more subtle, since the R dependence in both the soft 
and jet functions is relevant at leading power in r and since the dependence on R in both 
functions is singular as R —> 0. To remove it, we need a subtraction. To construct the 
subtraction, first recall that the general amplitude-level factorization proof in [FS2] applies 
to any operator, including one composed of Wilson lines. In particular, collinear factorization 
for a Wilson-line operator implies 


s = S R ® J* ® (34) 

where S = S R=oc has no angular or energy restriction and the eikonal jet function is 
defined as 

J*i W = £ I dn x 

X ! J 

The eikonal jet function differs from the jet function in Eq. (26) in that V/ replaces the field 
-0. Note that the measurement function in the eikonal jet function is the power-expanded 
version, 5{p — p + ), rather than 0 T (p). This is consistent with Eq. (34) since the phase space 
in the eikonal jet function is restricted to be in a cone. 

Explicitly, to order a s , we find 


(WI y/Wi |o) 


(0| Y^W\ |0) 


Kt-Px^) 
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(36) 


Comparing Eq. (36) to Eq. (28), we see that the cu and R dependence in J ei £ is the same as 
that in J Rj . This is expected, since the only IR-sensitive difference between the two is in 
the collinear-but-not-soft region of {X\ Y{W\ |0) and {X\ iJ;W\ |0). In this region, there is a 
complete cancellation of real and virtual graph for both functions, hence both are IR-hnite. 
Note also that there are no \ poles in the eikonal jet function. These double UV poles in 
the regular jet function come from virtual graphs. In the eikonal jet function, the virtual 
graphs in the numerator and denominator of Eq. (35) are identical and hence cancel in the 
ratio to order a s . The lack of \ poles also implies that there are no Sudakov double logs in 
the eikonal jet function. 
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Now, if we convolve both sides of Eq. (25) with the eikonal jet functions and use Eq. (31) 
and Eq. (34), we get 

^ ®J R2 ^HxS® J Rl ® J R2 (37) 

At this point, no object in this leading-power equivalence depends on A and the R dependence 
on both sides is only in the jet functions and eikonal jet functions. We still must have R 
small though, since there are power corrections in R on both sides. 

Finally, we want to remove the i?-dependence completely. Let us call a jet function 
with no restriction on R an inclusive jet function and denote it by .P. Removing the R 
introduces additional unphysical singularities collincar to the Wilson-line direction tj which 
are not regulated with the off-shellness regulator. We must introduced another regulator for 
these singularities, so we use the A-regulator [66], Sj, as shown in Eq. (51). To order a s we 
find for the inclusive jet function 


T ,, . . oi s Cf / n 2 \ £ f / 2 3 7 7r 2 


+ 5 (r) ( — In uj + 2 In uj In Sj + In 2 to j — f 2 In Sj + ^ 


' 1 ' 
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(38) 


Similarly, for the inclusive eikonal jet function we find 


4.m = *m + 2 £r($)'w 
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(39) 


Note that the Sj dependence associated with the Wilson line direction is identical in the two 
inclusive jet functions. 

Next, note that since Rj does not contain the jet direction the only leading-power con¬ 
tributions to the jet function from this region are soft. 1 Thus, we can apply the general 
amplitude-level factorization theorem to the operator Wj to get 


,P ^ 


J 


Rj 


tP 

Aik 


(40) 


This equation can be verified at 1-loop by comparing Eq. (38) with the combination of 
Eqs. (28), (36) and (39). Similarly, 


A = 4; ® 4k (4i) 

1 One might be concerned about collinear singularities associated with the Wilson-line direction tj. How¬ 

ever, since the measurement function forces p ■ n = t, at small r radiation cannot be collinear to both tj 1 
and the jet direction rP. Thus, radiation collinear to tj 1 cannot contribute at leading power. 
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( 42 ) 


Therefore, convolving both sides of Eq. (37) with and .J^ gives 

^ <8 J e \ k (8) J e 2 ik = H x S®J l ®J 2 

CLT 

In this final form, all the dependence on A or R has been explicitly removed. 

Finally, we want to isolate from Eq. (42). To do this, we use that convolutions map 
to products in Laplace space. Taking the Laplace transform, Eq. (42) translates to 


da 

dr—e 

dr 


— 1ST 


= H 


S{y)J l (y)J 2 (v) 


(43) 


This form is in agreement with previously known expressions in the literature [20,22], 


3.2 Jet broadening 


The above discussion shows how the phase-space cutoffs separating collinear and soft radia¬ 
tion as well as the UV phase-space cutoff can be removed in the factorization formula for a 
particular observable (thrust). The derivation easily generalizes to many other observables. 
The key general property that was used is that the vanishing limit of the observable forces 
the phase space into the IV-jet configuration at leading power. This allows the factorization 
theorem in Eq. (1) to be used to factorize the matrix-element squared in the full distribution. 
It also ensures that the dependence on the phase-space cutoffs is power suppressed, once the 
eikonal jet functions are included. For observables whose measurement function is not linear 
in each sector, the integrals will not be a simple convolution. 

For a marginally different example, consider jet broadening [46-50]. (Total) jet broaden¬ 
ing acting on a state \X) with particles of momenta p^ has the eigenvalue 


b(X) 



Mb(p) = \p± 


(44) 


where p± are the components of the 3-momenta of the particles transverse to the thrust axis. 

In the SCET literature, jet broadening is considered a SCETn observable because soft 
emissions which are hard enough to recoil against collinear emissions contribute to jet broad¬ 
ening at leading power, while they are subleading power for thrust. More explicitly, for thrust 
only the small component of momentum p + = n ■ p contributes (for particles going in the 
n hemisphere). Thus collinear momenta, with ( y p~ 1 p + 1 p ± ) ~ Q(1,A,A 2 ) and soft momenta 
with p ~ Q A 2 contribute at the same order. Soft momenta scaling like p ~ Q A give a 
power-suppressed contribution to thrust. For jet broadening, p± is measured. So collinear 
momenta contribute p± ~ Q A and therefore soft momenta scaling like p ~ Q A are relevant 
at leading power making jet broadening a SCETn observable. 

From the point of view of the factorization as set up in [FS1] and [FS2], the soft scaling 
is unrelated to the collinear scaling. That is, the amplitude-level factorization formula, 
Eq. (1) holds for any relationship between the soft-scaling parameter, A s , and the collinear- 
scaling parameter, A c . SCETi corresponds to \ s = A 2 and SCETn to A = A c . The relevant 
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implication of the soft and collinear momenta having p± components of the same order is 
that configurations where soft particles recoil against collinear particles must be accounted 
for in the factorization theorem. The result is that the factorization formula has the form 


% = H f db s db 1 db 2 d 2 pid 2 p2J 1 (b l ,Pi)J 2 (b 2 ,P2)S(b s , - p j 1 , -p 2 )S(b ^ b s -b x — b 2 ) (45) 

We can write this heuristically as 


da 

~db 


= H x J 1 ® J 2 ® S 


with the understanding that ® for jet broadening refers to the double convolution in Eq. 
With phase-space restrictions, these jet functions are given by 


./*■(!>,pP = £ 
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6lb-b(X 1 ))S(Q-p + Xi )S(p ± -pj Ci ) 


(46) 
(45). 

(47) 


and the soft function by 


s AR (b,Pi,Pi ) = 


£ / 


(A'.lpniO) sib - b(X,)) [S(pi - P^J + S(pi - f x J] 


(48) 


where p Xs is the net _L momenta in the left hemisphere and p 2 Xa is the net _L momenta 
in the right hemisphere. As with thrust, these phase-space restricted functions will have 
overlapping UV-IR divergences and unwieldy dependence on the cutoffs R and A. However, 
as with thrust, we can convolve both sides of Eq. (46) with eikonal jet functions to get a 
factorization formula with only objects with no phase space cutoffs. The result in Laplace 
space is 


db^e~ vb = H 
db 


'. 21-21 J{*L,v)J{xR,v)S{xi,xji,v) 

Q X Q X i 




Je ik(xf, v) J e ik(x^, v) 


(49) 


Here, v is the Laplace-conjugate variable to b and xj~ and x ^ are the Laplace conjugate 
variables to p ^ and p ^ respectively. 


3.3 Comparison to other approaches 

We have seen how phase-space cutoffs can be removed for certain inclusive observables if 
the double counting in the soft-collinear region is removed with an eikonal jet function. In 
this section, we would like to emphasize some conceptual differences in our derivation and 
previous ones and contrast with the literature. 

First of all, it is easy to compare our results to those in traditional QCD, where the eikonal 
jet function first appeared. The final factorization formulas are identical. One difference is 
that in the early literature the eikonal jet functions were subtracted from the soft function 
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rather than the jet function. From the point of view of the final formula, there is no difference. 
However, conceptually our analysis makes it clear that the eikonal jet function should be 
subtracted from the jet function rather than the soft function. Indeed, the soft function is 
by itself infrared finite while the naive inclusive jet function (without the subtraction) is not. 
As shown explicitly in Eq. (38), the infrared divergences do not cancel between the real and 
virtual graphs for the jet function. With the subtraction the jet function is a well-defined 
and infrared safe object. The fact that the subtraction is more naturally applied in the 
collinear sector was also appreciated in [22], 

The comparison to SCET is perhaps more illuminating than the comparison to traditional 
QCD. In the early days of SCET, calculations were mostly done in dimensional regularization 
(DR) and the overlapping of soft and collinear phase-space regions were not much discussed. 
In retrospect, it is easy to see why the correct answers result in DR without a subtraction: 
the eikonal jet functions, as in Eq. (35), give scaleless integrals in DR and thus formally 
vanish and can be ignored. 

It is natural to be somewhat uncomfortable with setting scaleless but IR and UV divergent 
integrals to zero in DR. The mathematical justification notwithstanding, it is dangerous from 
a practical point of view if one hopes to extract an anomalous dimension from the poles in 
the jet and soft functions at d — 4. The only way it will work is if the object one computes 
is infrared finite. For infrared hnite objects, all the poles are by definition UV. So setting 
—- — = 0 has no effect. As we have shown, the subtracted inclusive jet function is IR 

^UV £IR 

hnite, so practically, one can ignore the subtraction in DR. Morally, though, to do this one 
must be able to show that the jet function is IR hnite. Without the subtraction it is not. In 
this respect, the success of SCET in pure DR was somewhat accidental. 

The missing subtractions were understood in the classic paper on zero-bin subtraction by 
Manohar and Stewart [21], These authors showed the the proper derivation of the effective 
Lagrangian for SCET involves binning the momenta into collinear momenta in different 
directions and soft momenta. The zero bin in each collinear sector should be formally 
excluded. In [21] it was shown that this exclusion amounts to the subtraction diagram-by¬ 
diagram of the soft-limit of the collinear momenta. In [22-24], this subtraction procedure 
was shown to be equivalent to the eikonal jet function subtraction of traditional QCD. 

A somewhat different perspective comes from the method-of-regions multipolc-expansion 
approach to SCET [25-27], as recently reviewed in [28]. In this approach, the soft-collincar 
subtraction and the extension of the soft integrals to A = oo are not discussed or needed. 
The basis of the argument is that the non-analytic dependence on the observable in each 
region is independent of the possible phase-space cutoffs. Since all the physics is in this non- 
analytic dependence, one can remove the cutoffs without consequence. For more details, 
see [25,28], 
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4 Conclusion 


In this paper we have presented two new results. First, we have given a recursive formula for 
constructing amplitudes which agree with full QCD at leading power to all-loop order in any 
soft and collinear limit. The subtracted amplitudes we describe are matrix elements of fields 
and Wilson lines. Unlike with the amplitude-level factorization formula in [FS2], one does 
not have to specify whether the particles are soft or collinear ahead of time: the subtracted 
matrix elements will be correct in any limit. Although the amplitudes appear simpler than in 
full QCD (for example, the only interference effects from different directions involve gluons 
emitted off Wilson lines), it remains to be seen whether they can be integrated simply to 
provide a productive subtraction scheme. In our derivation of this formula, extensive use 
was made of the proof of factorization in [FS2], 

Second, we showed how phase-space cutoffs can be removed when integrating a factor¬ 
ized amplitude squared against the measurement function for certain inclusive observables. 
Removing the cutoffs does two things: it overcounts the soft-collinear region and adds UV 
divergences to the phase-space integrals. These two effects can be compensated for by inte¬ 
grating the full QCD distribution against an eikonal jet function. This convolution can be 
easily disentangled, at least for thrust, jet broadening and angularities. This extends the 
results of amplitude-level factorization from [FS1] and [FS2] to the level of observables. 

In our presentation, we have included explicit 1-loop expressions for soft and jet functions 
with cutoffs and for the eikonal jet function in a regularization scheme which separates UV 
from IR. These expressions confirm generally the qualitative analyses that we have presented 
of the UV and IR structure of the integrals. 

Although our final factorization formulas are not new, we believe our derivation is sys¬ 
tematic and rigorous. We hope that the step-by-step procedure we have presented will be 
useful in future studies of factorization, where subtleties abound. 
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6 Appendix 

In this appendix, we report the explicit computations of many of the cut-off dependent 
objects that appear in Section 3.1. All expressions are computed at one-loop order in QCD 
for the observable thrust. 
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6.1 Regularization schemes 


We will use dimensional regularization to control the UV divergences: we analytically con¬ 
tinue to d = 4 — 2 e dimensions, with e > 0. 

To regulate the IR divergences, we take the outgoing external fermion lines to have an 
offshellness of Q 2 uj. Consequently, the propagator in Y) will look like 


n- 


u, 

Pj 


n 3 ■ k 


-)•-(50) 

Pj ■ k + Wf- 

It is helpful when using an offshellness also to slightly modify the measurement function for 
thrust, so that virtual graphs still contribute only at r = 0. We can do this by replacing 
S(r- ^ p~ ) ->• 5{t - ±p~ + uj). 

In scaleless integrals involving Wilson lines, the offshellness may not completely control 
all the IR divergences. Thus, in addition we use the A regulator [66] for the collincar Wilson 
line, so that the propagator in the Wj will be shifted by -r^— . More precisely, we define the 

J Zj -pj 

dimensionless parameter 5j = A /(tj ■ pj) and shift the eikonal propagators as 

11 1 


tj-k 




ij-k+tpt~ ij ■ k + 8j(tj ■ pj) 


(51) 


With these two IR regulators, all the 1/e poles in the following expressions correspond to 
UV divergences only. 

For simplicity, in the following sections we choose the direction of the collinear Wilson 
line, W\ to be t 1 / = n/, = (1, —n) and that of Wo to be t% = n 1 / = (1, n), where n is the thrust 
axis. Then, in light-cone coordinates 

k M = ^k~rii + ^jk + ri2 + kj_ = {k ~, k + , k±) (52) 

While it is not impossible to do the integrals for generic choices of t/ and t%, the phase space 
integrals inside the cone become significantly more complicated. 


6.2 Jet Functions 


We will denote a jet function restricted to a cone of size R by J R . In our notation, R is a 
rapidity-type cone angle: a particle is in the n-cone if k + < Rk~. R is related to the cone 
opening angle 9 by R = tan 2 1. 

At order a s , the virtual contributions to the jet functions are of course independent of 
phase space restrictions and therefore the same for the restricted or unrestricted jet functions. 
The purely virtual contributions come from evaluating 
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where the sum over spins is implicit. 

There are 3 virtual graphs contributing to (pj\ i^Wj |0) at 1-loop. The self-energy graph 
on the Wilson line exactly vanishes; the vertex correction and quark self-energy graph sum 
to 
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In the last step, we expanded in Sj and to and dropped the 0(Sj) and 0(uj) terms. 
To one-loop, the denominator factor in the jet function evaluates to 
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We note here that the Sj dependence is identical in Eqs. (55) and (56). Thus virtual IR diver¬ 
gences introduced by the Wj, which are regulated with the A-regulator, cancel in [J(r)] V i rt . 
to one loop. Explicitly, 
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Now let’s look at the real-emission diagrams. With emissions restricted to a cone of size 
R , the real-emission contributions to the jet function at order a s are 
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d p + f dp 


real 


27T / 27T 


{2tt)5(p - Q ) 


d d k 1 / 2 [p — k ) d — 2 p + — k + 
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p-i 


2nS(k + k- - kl)d(k + )6(k- 


J (27T ) d p + \k + +p~6 1 
x 2nS(^(p + — k + )(p~ — k~) — k\ — Q 2 u?j9(p + — k + )9(p~ — k~)S(r — p + /Q + c o) 
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It is helpful to write the result as 


J Rj (t) 


oc,C 


s'-F 


-real 2ll T(l — e) \ Q 2 J T + UJ T e 


Anjl 


2 \ £ 


X R Hr) 


(59) 


and we find 

z R n 


t)= dx 
Jr/R 


2(1 -x) 
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(60) 


where Fi(a, /3, (3'] 7 ; x, y) is the Appell hypergeometric function and B(z]a,b) is the Incom¬ 
plete Beta Function. These real emission graphs are of course UV finite, so we can simply 
set 8 = 0 . 

To expand the result for small c 0 , we can use that in the limit uj —> 0 


r + u 


X(t) = S(t) 


dr' 


Jo t' + u 


X(t') 


+ 


X(r) 


(61) 


so that 
J r Jt) 


real 


a s Cp 

2 tt 


5(r) 


- ln 00 — 2 ln uj ln R + In 2 uo + 0(R ) 


+ 


R 2 


(62) 


Eq. (62) shows that the jet function in cone does not have log^- singularity. Note, however, 
that this jet function does have double logs of u, coming from soft-collinear region; these 
will cancel against the virtual soft-collinear singularity in Eq. (56). The R dependent terms 
logo;log R come from soft emissions at the cone edge. These do not cancel against any other 
contribution to J R , but will cancel against either S R of the eikonal jet function Jj k . 

Next, consider the inclusive jet function, JX This is defined identically to J R , but 
without the phase space restriction. The virtual contributions are the same. The real 
emission contributions at order a s are the same without the phase space restriction on the 
last line of Eq. (58). Using the same notation as in Eq. (59), we find 


V = / dx 


2(1 -x) 
x + Sj 


+ (1 - e)x 


x(l — x) 


4- 1 +£ v / vrT(l - e) 

^■r(l-e) 


— (3 + 8)Sj +4(1 + SJ) 2 Fi( 1 , 1 — 8,2 — 28, ——) 


- -- + 2(1 + 6,) ln(l + —) 


(63) 

(64) 
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We have set e = 0 on the last line since the real emission contribution to the inclusive jet 
function, like with J R , is UV finite. Using Eq. (61) we then find 


J j (r) 


real 


| <5(r) hi cu + 2 In u In Sj ) + 


— ( —2 In 5; 



(65) 


The phase space integral in the inclusive jet function contains single log of u>, indicating pure 
collinear singularity. It also contains a double IR singularity, log u log Sj, coming from the 
soft-collinear region. 

Next, we consider the eikonal jet function, with the held 0 replaced by a soft Wilson line. 
The virtual contributions are given by matrix elements similar to Eq. (53): 

2 


Je ik(^") 


dp H 


virt. 


2 vr 


( 0 | Wi| 0 > 


( 0 | y\W\ | 0 ) 


27rh (p + - Qu?j 9(p + )S(r + u) 


= 


( 66 ) 

(67) 


Since there is an exact cancellation between numerator and denominator, the purely virtual 
contribution is S(r) to all orders. 

Using the same notation as above, for the phase-space restricted eikonal jet function, we 
End 



S~‘B 



?-21 n (1 + 0 + 0 ( £ ) 


( 68 ) 


Note that now there is a UV divergence from the integral up to infinite energy in the k 
direction. This turns into 
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’2vr 2 2 , ' 


ri (2 \] 
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For the unrestricted eikonal jet function 
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^•eik ” 


r»oo r\ 

dx 


x + Si 


= S-‘ 2r(e)r(l-e) 


(69) 

(70) 


and 


4k ( r ) 


_ a s C F X £ 

real 27T V^ 2 . 

x i <5( r ) f~^ -In o; + 2 In uj In Sj + In 2 to ) + 


— (-2 In 5.; — 2 In r 


T \ £ 


(71) 


Finally, the jet functions with real and virtual contributions combined are given in 
Eqs. (28), (36), (38) and (39) for J R \ ,J ^ , J j and J 3 eik respectively. The structure of 
IR singularities in virtual and real contributions to the various objects are listed in Table 1. 
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Tabic 1: IR singularities in virtual and real phase-space integrals at order a s . Note that the 
jet functions and eikonal jet functions are all IR divergent but the IR divergences cancels in 
their difference. Indeed, the subtracted jet function is infrared safe and has no dependence 
on ca or Sj. Note that logo;log R and logo;log Sj terms describe soft-collinear singularities. 
After subtracting the eikonal jet function from the inclusive jet function, these singularities 
drop out. 
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6.3 Soft functions 


Here we focus on the soft function for thrust, defined as 


S(r) = £/ 


dU Xs 



( 72 ) 


with Q r (p) the measurement function for thrust defined in Eq. (22) and X s the appropriately 
phase-space restricted set of states. 

The virtual contribution to the soft function comes from the matrix element (0| Y^Y 2 |0) 
and is independent of phase-space restrictions. To order a s , 


(0| Y'Y 2 |0) = 1 + g 2 s p 2e C F 
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(73) 


(74) 

(75) 


We see that the virtual contribution has double poles in the UV and IR. In pure dimensional 
regularization, these double poles would exactly cancel and the virtual contribution exactly 
vanishes. 

Next, we compute the real-emission contributions. First, we consider the inclusive soft 
function. At order a s , the real emission contribution is given 


h(T)],„, = a g yc F 


d d k 


(2n) d \k + + Quj k~ + Qlu 


2t i5{k + k- - kl)e(k + )6(k ~) 
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6(k - k + )5(r - ^—) + 6(k + - k )<5(r - ^—) 


As with the jet function, we write the result as 
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a K C, 


s'-'F 
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(76) 


(77) 


We will use this parameterization for the inclusive and phase-space restricted soft functions. 
For the inclusive soft function, we find 


XAt) = 


dx 


so that, using Eq. (61), 
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real 
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adding the virtual contribution, we get to order a, 


S(r) = S(t) + 


OI s Cf 

7T 





( 80 ) 


which is IR finite (no dependence on uj). The UV poles are removed by renormalization. 

The inclusive soft function is IR finite, since all uj divergences cancel between real and 
virtual diagrams, as we can see in Eq. (79) and Eq. (74). As with the jet functions, 0R es , 
represents the phase-space restriction. For the soft function in the ball, the phase space 
restriction is ©a = 9( A — k~)9(k~ — k + ) + 9( A — k + )9(k + — k~)\ for the soft function outside 
the cones, ©^ = 9(R — k~ / k + )9(k~ / k + — 1 /R); finally, for the soft function in the ball and 
outside the cones, © A jj = ©a©^. 

For the soft function phase-space restricted to radiation outside of cones around the jet 
axes (but with no cutoff A on energy), we find 

= 2 In + uj'j — 2 In (r + lS) + 0(e) (81) 

Note that this soft function is not UV divergent, since at finite r there is an implicit energy 
cutoff E < |(1 + jj), as discussed below Eq. (32). It leads to 
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71 
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5(t) ( — — + 2 In uj In R + 0(R,uj) ) + 


-Ini? 

T 


(82) 


Note that the IR-divergent In a; Ini? term does not cancel against the virtual correction. 
Thus the soft function outside the cones is not infrared safe. 

For the soft function phase-space restricted to a ball with energy less than A, but no 
angular restriction, we find 



2 ln(r + uj) 


(83) 


This is also UV finite, as expected. To get the soft function outside the cones with energy 
less then A, we can simply combine these two expressions with ^-functions: 


I» AR (r) = - a) +I, r 9(a - jj) 


(84) 
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which leads to 



5(t)+C f — 
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— — + 2 In uj In R + O(R) 
o 


-Ini? 

T 



a s 


-C F — 


7r 
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(85) 


Combining this with the virtual contribution leads to Eq. (29). Note that for r sufficiently 
small (r < A R), X AR and hence S AR does not depend on A. In particular, at leading power 
in r we have 

S AR ~ S {A=oo)R (86) 

as in Eq. (31). 

To repeat, we fold that the inclusive soft function is IR finite while S R and S AR are not, 
due to an incomplete cancellation of IR singularities between virtual and real contributions. 
Indeed, S R contains a In uj In R term, coming from the real-soft region, and a In 2 a; term, 
coming from the virtual-soft-collinear region. These singular terms exactly match with those 
in J^, according to Eq. (82) and Eq. (69). 
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